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SECTION  I 


INTRODUCTION 

V 

Two  separate  but  related  efforts  are  described  in  this 
report.  The  first  and  more  significant  effort  is  a demonstration 
of  the  capability  of  the  patch  zoning  method  for  numerically 
solving  the  magnetic  field  integral  equation  (MFIE) . This  is 
accomplished  by  comparing  experimental  data  with  numerical 
computations  of  the  current  density  induced  on  metallic  struc- 
tures. The  most  immediate  impact  of  this  capability  is  its 
effect  on  determining  the  potential  and  limitations  of  experi- 
mental procedures. 

The  structures  that  were  treated  numerically  and  experi- 
mentally are  an  aircraft  model  in  free  space  and  a metallic 
cylinder  of  finite  length  in  free  space  and  also  above  and 
parallel  to  a metallic  ground  plane.  The  experiments  were 
performed  at  the  University  of  Michigan  by  Valdis  Liepa  and 
his  associates  under  contract  to  the  Air  Force  Weapons  Labora- 
tory. 


The  secondary  effort  is  a derivation  of  the  electric  field 
integrodif ferential  equation  (EFIDE)  and  the  MFIE  for  electro- 
magnetic pulse  external  interaction  with  perfectly  conducting 
bodies  above  a finitely  conducting  half-space.  The  relation- 
ship between  the  two  efforts  is  that  the  equations  are  repre- 
sented in  such  a manner  that  the  computational  results,  which 
are  compared  with  experimental  data,  are  readily  seen  to 
correspond  to  the  numerical  treatment  of  special  cases  of  the 
MFIE  finitely  conducting  ground  equations. 

The  derivation  of  the  EFIDE  and  the  MFIE,  which  include 
lossy  earth  effects,  utilizes  the  explicit  representations  of 


Green's  dyadics  given  by  Tai  (ref.  1)  after  appropriately 
accounting  for  missing  terms  (ref.  2) . The  results  obtained 
by  the  straightforward  use  of  Tai's  dyadic  representations 
contain  infinite  sums  of  infinite  integrals.  We  introduce  a 
procedure  that  primarily  consists  of  a coordinate  change  that 
reduces  these  sums  to  finite  sums  of  well  studied  Sommerfeld 
integrals.  The  reduction  of  the  infinite  sums  to  finite  sums 
was  not  unexpected  since  we  could  have  obtained  the  same 
results  by  appropriately  combining  vertical  and  horizontal 
dipole  solutions. 

To  enhance  the  utility  of  our  equations,  we  change  Tai's 

f\j 

notation  to  that  of  Banos  (ref.  3)  after  the  infinite  sum 
reductions.  We  do  this  for  two  reasons.  First,  Banos  presents 
a detailed  theoretical  investigation  of  the  resulting  Sommer- 
feld integrals.  The  other  reason  is  that  the  work  of  Lytle 
and  Lager  (refs.  4,  5)  uses  Banos  notation  and  it  contains  com- 
puter codes  for  the  evaluation  of  the  Sommerfeld  integrals. 

The  recent  work  of  Haddad  and  Chang  (ref.  6)  also  contains  both 
theoretical  and  numerical  work  related  to  the  evaluation  of 
the  Sommerfeld  integrals. 
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SECTION  II 


PRESENTATION  OF  UNSIMPLIFIED  REPRESENTATIONS 
FOR  E(r)  AND  H(r) 

For  the  purpose  of  this  analysis  it  is  convenient  to  divide 
our  volume  of  interest  into  two  regions  (fig.  1),  Vp  and  Vc. 

Vp  is  the  volume  bounded  by  the  surface  of  the  object,  S,  the 
interface  surface,  Sg,  and  the  upper  hemisphere  at  infinity, 
SFjo.  Vp  includes  the  volume  Vj  over  which  the  source  J(r)  is 
defined.  is  the  semi- inf inite  volume  bounded  by  Sg  and  the 


lower  hemisphere  at  infinity 
region  are 

The  equations 

satisfied 

in  each 

V x E,  (r ) = io)u 
—I  — o 

Mi  (r) 

r e VF 

(1) 

7 x H,  (r)  = - itoe 
—I  — i 

3 ^(r)  + J(r) 

I € VF 

(2) 

and  J(r)  = 0 unless  r e Vj . 

The  equations  in 

Vc  are 

7 x Eg(r)  = icowo 

M2  <r-> 

r € Vc 

(3) 

V x H?  (r)  = -iuje 

E2(r) 

r € Vc 

(4) 

where  e is  a complex  function  of  cj  [i.e.,  e = eR(cu)  + ie^Co))]  to 
account  for  the  fact  that  the  half  space  is  lossy  and 
frequency  dispersive.  Combining  equations  1 and  2 we  obtain 

Lp  E^Cr)  = iu)uQ  J(r)  (5) 

where 

L = y x v x - k2  (6) 

r O 

7 

t 


Figure  1.  Geometry  for  Electromagnetic  Scattering 
from  a Perfectly  Conducting  Body  Situated  Above 
a Finitely  Conducting  Half  Space. 
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V?  ’ 


and 


Similarly  combining  equations  3 and  4 we  obtain 

Lc  E2(r)  = 0 (7) 

where 

Lc  = 7 x V x -k2  (8) 

and 

k2  = w2uQe.  (9) 

Next  we  introduce  the  Green's  dyadics  that  satisfy  the  equations 
LF  Si  (£.£<,)  •!«(£.  rc)  £.  £q  « VF  (10) 

and 

Lc  G2(r,r  ) ■ T 6 (r  - r ) - f VC 

C 2 ^ ^ lo  € VF-  Cll) 

Next  we  use  the  dyadic  identity  for  an  arbitrary  vector  a and 
dyadic  A found  in  reference  7.  That  is 
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*•**■  ■***-.■  M -f 


(V  x a))  • A - a 


(V  x 


dV 


(n  x a)  • (V  x A)  + 


(n  x (v  x a))  • A f dS  (12) 


where  V is  an  arbitrary  volume  bounded  by  the  closed  surface  S 
which  has  an  outward  normal  n.  By  adding  and  subtracting  either 

9 =*  9 = 

-k~a  • A or  -k  a • A the  identity  (eq.  12)  can  be  expressed  as 
either  of  the  following  two  identities. 

A - a • (LpA)J’  dV 

x a)  • (V  x a)  + (n  x (V  x a))  • A J-  dS  (13) 


acA)|  dv 

x A)  + (n  x (7  x 


a))  • A } dS.  (14) 


We  now  employ  equation  13  where  a = (r ) and  A = G^(r,r^)  and 
V = . Using  equations  1,  5 and  10  we  obtain 


S<Eo> 


+ -S  + — SB1 


— SF° 


(15) 


10 


*— 


where 


= / (ia)Wo)(j(r)  • 31(rfr<J))  dv  (16) 

VT 


ls  = iwuo  / (ng  (r)  x H(r))  • G-^r.r^  dS  (17) 

s 

ISFoo  = - /*  {(«-<£>  x Ii<£))  * (v  x ^1  (£•£<,)) 

*SF°° 

+ iwu0(nm(r)  x Hx  (r ) ) • ^(r.r^  dS  (18) 

I-.Bl  ‘ -/ * M£>)  • <’  » E1<E’Eo>'> 

SB 

+ iuu0(nB1(r)  x Hx(r))  • 51(r,r<))|dS  (19) 


and  n (r)  is  the  outward  normal  to  S.  Similarly,  employing 
s — __  __ 

equation  14  for  a = E2(r),  A = G2  (r,^)  , V = Vc  and  using 
equations  3,  7 and  11  we  obtain 


/( E2(r>  • (i  4(r  - rj) 


dV  = ISB2  + isc= 


V, 


where 


— SB2 


•/{( "B2(^  * -2(^  * (V  * •!<>>) 


+ i^o(nB2  x * ^(r.r^)}  dS 


(20) 


(21) 


t 


11 


r 


Iso  * ■/  {("-  * • (v  x 


sc« 


+ i“>U0  (”«  * —2  dS' 


(22) 


we  note  that  ICP 
— SF°° 

—SC  00 


We  now  simplify  equations  15  and  20  by  employing  the  boundary 
conditions  satisfied  by  the  fields  and  the  Green's  dyadics . First 

0 as  a result  of  the  radiation  condition  and 
0 as  a result  of  the  exponential  decay  due  to  losses  (the 
radiation  condition  would  be  sufficient  if  e were  purely  real). 
Next  we  note  that  the  volume  integral  in  equation  20  is  zero 
because  the  integration  is  over  and  according  to  equation  11, 
r^  € Vp.  Combining  these  results  we  can  write  equations  15  and 
20  as 


E(r) 


+ icou 


A 


(r) 


Gl(£.Eo> 


dS  + — SB1 


(23) 


We  now  return  to  equation  17  which  combined  with  the  information 
contained  in  equation  26  enables  the  interpretation  (and  evalu- 
ation for  plane  wave  excitation)  of  E^Cr^)  without  explicitly 
performing  the  indicated  integral  over  Vj.  From  equation  16  we 
see  that  E^tr^  is  independent  of  the  surface  S.  From  equation 
26  we  see  that  by  imagining  S to  vanish,  E^r^  can  be  inter- 
preted as  the  total  field,  incident  plus  scattered,  due  to  the 
lossy  half-space  with  S absent.  For  an  incident  plane  wave,  ET ( ) 
would  be  simply  expressed  (algebraically)  in  terms  of  the  Fresnel 
reflection  coefficients. 

We  now  rewrite  equation  26  with  the  following  change  variables 
to  conform  to  standard  notation,  i.e.,  r^  is  now  denoted  as  r and 
r is  now  denoted  as  r ' . The  resulting  equation  is 

dS'  (27) 

which  can  be  rewritten  as 

E(r)  = ET(r)  + iouQ  j* ^ (r ' , r)  • Js(r')  dS 1 (28) 

s 

where  'v  denotes  transpose.  We  now  use  the  fact  that 

^(r’  ,r)  = 5^U)  (r.r*  ) (29) 

where  (r,r’ ) is  the  notation  used  by  Tai  (ref.  1)  and  has 

exactly  the  same  meaning  implied  in  that  book.  Combining  equations 
28  and  29  we  have 

E (r)  - ET(r)  + iwu0  f C<n)(r,r')  • Jg(r')  dS'  . (30) 


E(r)  = 


ET(r) 


1U)U 


J (r' ) 
— s — 


G1(r' , r) 


; 


13 


Up  to  this  point,  all  of  the  work  presented  has  been  tutorial; 
now  we  make  our  first  significant  point.  As  pointed  out  by  Tai 
(ref.  2),  the  dyadics  presented  in  his  book  lack  necessary  terms. 
His  method  of  correcting  these  dyadics  leads  to  expressions  that 
do  not  explicitly  exhibit  the  best  form  for  subsequent  numerical 
treatment.  For  the  lossy  half-space  dyadic  (and  others  as  well) 
we  can  express  the  dyadic  as  the  sum  of  two  terms 

£3(11)(r,r')  - S0(r,r')  + 5^1)(r,r')  (31) 

and  the  omitted  terms  were  omitted  only  from  G (r,r‘),  the  free 
space  Green's  dyadic.  By  correcting  Go(r,r' ) by  supplying  the 
missing  L related  dyadic  terms  (or  equivalently  the  missing  6- 
function  term)  we  would  end  up  with  a representation  for  Gq( r,r’) 
that  would  not  be  as  useful  for  numerical  purposes  as  the  standard 
representation  that  does  not  contain  the  L,  M,  and  N related 
dyadics.  The  following  two  observations  are  the  basis  of  this 
claim : 

1)  The  expanded  version  of  Go(r,r’)  will  always  contain  two 
different  representations  that  require  separate  numerical  treat- 
ment, depending  on  whether  one  of  the  spatial  observation  point 
coordinates  is  larger  or  smaller  than  the  corresponding  integration 
point  coordinate. 

2)  The  derivation  of  the  integral  equation  in  the  magnetic 
field  case,  or  the  integrodif f erential  equation  in  the  electric 
field  case,  requires  that  the  limit  be  properly  treated  as  the 
observation  point  approaches  the  integration  point  on  the  surface 
S.  By  expressing  Go(r,r‘)  in  the  standard  closed  form,  we  can 
make  use  of  existing  analysis  to  treat  this  limit.  With  this  as 
background  we  have  the  following  representations: 

G0(£,r’)  - (l  + w)g  (32a) 

k 
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where 


g 


exp  (ikQ|r 


(32 


and 


G<U> 

3s 


(r.r') 


00 


Oo 


ME  E 

o n*0  a=e , o 


(2 


5 )iaM  ,(h,)M'  ,(h,) 
no  ( — anA  1 —anA  1 


where 


+ “WVSinxOV} 

«enx<hl)  * 5 * [V*r)stn  "♦  cihlZ  *] 
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hl  ' h2 
hi  + h2 


b = 
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k2hl  ~ koh2 

k2hl  + koh2 


2 

to  u e 

o o 


(3 


(3 

(3 


(2 

(2 


k?  = a)2u  e;  e is  the  complex  permittivity 
^ ° of  the  lossy  half-space  (' 


and  6no  is  the  Kronecker  delta  function.  Implied,  by  the  explicit 
relationships  presented  in  equation  33,  is  a coordinate  system 
having  its  origin  at  the  interface  between  the  serai- inf inite 
lossless  half-space  and  the  serai-infinite  lossy  half-space.  The 
usual  cylindrical  coordinate  system  is  employed  with  z being 
measured  as  positive  as  the  distance  from  the  interface  is  increased 
into  the  lossless  medium.  In  summary,  equations  28  through  33 
present  an  explicit  representation  for  the  electric  field,  E(r) , 
off  the  surface  S once  Jg  (r ) has  been  determined.  A discussion 
concerning  the  determination  of  an  appropriate  equation  having 
Jg (r ) as  its  solution  will  be  given  in  subsequent  sections. 

The  equations  that  present  an  explicit  representation  for  H(r)  off 
the  surface  S,  once  Jg  (r)  is  determined,  are  arrived  at  by  taking 
the  curl  of  equation  28  and  employing  equation  1.  The  resulting 
representation  is 

')  d S'  (34) 

S 


H (r ) = HT(r) 


j K(r,r’ 


J 

— s 


(r 


where 


K(r , r ' ) = V x Gj (11) (r,r' ) (35) 

and  employing  the  same  arguments  that  led  to  the  interpretation 
of  Ej(r)  without  the  explicit  evaluation  of  equation  16,  we  can 
interpret 


Sr®  * lir  v x &r<£>  <36> 

o 

as  the  total  magnetic  field  ( (r ) ) , incident  plus  scattered,  due 
to  the  lossy  half  space  with  S absent.  For  an  incident  plane 
wave,  H.p(r)  would  simply  be  algebraically  expressed  in  terms  of 
Fresnel  reflection  coefficients. 

We  now  simplify  equation  35,  to  a limited  extent,  by  straight- 
forward substitution.  In  the  next  section  we  will  simplify  it 


a great  deal  more  and  discuss  the  numerical  benefits  of  the  form 
we  ultimately  obtain.  First  we  write 


K(r.r')  = Ko(r,r’)  + Kg(r,r') 

where 

= 7 x Go(r,r')  = Vg  x | 

and 

Ks  = V x G^l)  (r  , r ' ) 

with  S^^^Cr.r')  given  by  equation  33.  Equation  39  can  be 
fied  by  employing  equation  33c  and  the  relationship 


7 * N|n X Chl > * ko^nX<hl> 


The  expression  obtained  is 


ik 

Ks^'>  = 4 IT 


OO  CO 


E <2  - Sno>{ 


aN  (h, )M ' , (h  ' 
-an  X 1 -an  X 1 


n=0  a=e,o 


b— an  A 1 ^ —an  X ^ 1 T 


The  resulting  expression  for  H(r)  is 


H(r) 


HT(r)  + 


x Js(r’)(  dS’ 


+ 


Js(r’))  dS’  . 


(37) 

(38) 

(39) 
simpli- 

(40) 


(41) 


(42) 


17 


SECTION  III 


PRESENTATION  OF  THE  ELECTRIC  FIELD  INTEGRODIFFERENTIAL  EQUATION 
(EFIDE)  AND  THE  MAGNETIC  FIELD  INTEGRAL  EQUATION  (MFIE) 

IN  UNSIMPLIFIED  FORM 


The  representation  for  the  electric  field  off  the  surface  is 
given  by  combining  equations  30,  31,  and  32  to  obtain 


E(r)  - Ex(r)  + iu„0  J/s|IE.  - r'  |)  J.s(r')  dS' 


?/’ 

O S 


/' 


+ I Wg  • Jg (r ' ) dS’  + | e^1)(r,r')  • Js(r')  dS'j 


(43) 


We  take  the  cross  product  of  both  sides  of  equation  43  with  n(r) 
where  r is  the  point  on  S approached  by  the  r in  'equation  43 
which  was  a point  off  the  surface  in  that  equation.  We  also 
consider  the  behavior  of  the  second  integral  in  equation  43  as 
the  volume  observation  point  approaches  the  surface.  (The  behavior 
of  the  other  two  integrals  requires  no  special  treatment.)  The 
behavior  of  this  limit  was  treated  in  reference  8 and  we  can 
directly  use  those  results  by  employing  representation  31. 

Finally,  by  following  the  described  steps  and  employing  the 
boundary  condition 

n(r)  x E (r)  =0  r * S (44) 


we  obtain  our  unsimplif ied  EFIDE 
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X 


-n(r)  * (r) 


iwuo  n (r) 


!/ g (|r  - r’  |)  J8(r')  dS' 


LJs(r')  + Jc^^Cr.r’)  • Jg (r ’ ) dS’j 


(45) 


where  equivalent  representations  of  L are 


k^  LJg (r ' ) - 77  • Jg  (|r  - r'  |)  Js(r')  dS’ 

(46a) 

s 

- g (Ir  - r'  I)  V • J8(r')  dS  ’ 

(46b) 

- J^Vg  (| r - r'  |)  V’  • Jg(r')  dS’ 

(46c) 

- Jgll  - r*  | V*V  • Js(r')  dS* 

(46d) 

- vJ*Vg  (|r  - r’  I)  • J g (r ' ) dS’ 

(46e) 

and  the  principal  value  sign  is  only  employed  when  it  is 
required. 

The  MFIE  is  obtained  by  taking  the  cross  of  both  sides  of 
equation  42  with  the  n(r)  just  described.  We  also  consider 
the  behavior  of  the  first  integral  in  equation  42  as  the  volume 
observation  point  approaches  the  surface.  This  limit  is 
precisely  the  one  treated  in  the  free  space  MFIE  and  its 
result  is  known.  The  second  integral  in  equation  42  requires 
no  special  treatment  as  the  volume  observation  point  approaches- 
the  surface  observation  point.  Again,  this  is  a major 
benefit  of  the  Green's  dyadic  representation  31.  The 


■>:Y‘ 
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MFIE  resulting  from  taking  the  described  cross  product  and 
limit  is 


f(fi)  J (r)  = n (r)  x ^(r)  + J {n(r)  x (Vg  x J(r'))}  dS’ 


+ J{n(r)  x Ks(r,r’)  • J(r')}  dS ' 
S 


where  we  have  used  the  definition  (25)  and  dropped  the  sub- 
scripts on  J (r) . The  f (ft)  comes  from  the  limiting  process 
associated  with  the  first  integral  in  equation  42  and  is 


f(ft)  = 1 - ft/4ir 


where  ft  is  the  solid  angle  subtended  by  the  surface  S at  r. 
If  we  don't  choose  r at  a discontinuity  in  curvature,  then 
ft  * 2 -rr  and  f(ft)  assumes  the  value  1/2  which  is  usually  seen 
in  the  magnetic  field  integral  equation. 


SECTION  IV 


DERIVATION  OF  THE  TRACTABLE  FORMS  FOR  THE  EFIDE 

AND  THE  MFIE 

An  examination  of  the  EFIDE  presented  in  equation  45  and 
the  MFIE  in  equation  47  reveals  that  the  subsequent  numerical 
treatment  of  either  of  these  equations  requires  the  treatment 
of  an  infinite  sum  of  infinite  integrals.  This  is  the  case 
because  equation  45  contains  (£ ,r ')  defined  by  equation  33 

and  equation  47  contains  K (r,r')  defined  by  equation  41.  The 
objective  of  the  work  presented  in  this  section  is  to  derive 
new  exact  representations  for  G^^(r.r’)  and  K (r,r')  as  finite 
sums  of  demonstrably  convergent  integrals.  These  new  represen- 
tations when  substituted  into  equations  45  and  47  constitute 
the  tractable  forms  for  the  EFIDE  and  the  MFIE. 

The  method  of  obtaining  the  simplified  forms  relies  on 
the  interpretation  of  the  r and  r'  that  appear  in  G^^Cr.r') 
presented  for  the  first  time  in  equation  29.  After  the  change 
in  variables  and  transposition,  then  G^^(r,r')  satisfies  the 
equation 


V x v x G^11)(r,r')  - k°  G^11)(r,r’)  = IS(r  - r ' ) (49) 

where  the  derivatives  associated  with  the  curl  operator  are 
with  respect  to  the  r variation,  and  the  boundary  conditions  at 
infinity  and  at  the  interface  are  applied  to  r approaching 
infinity  and  r on  the  interface.  This  interpretation  means 
that  r'  is  not  required  to  vary  for  representations  33  and 
41  to  be  valid.  Even  though  r'  is  not  required  to  be  a 
variable  for  those  representations,  the  application  of  those 
representations,  for  example  equations  45  and  47,  do  require 
that  r'  be  a variable  point  that  has  no  particular  restrictions. 
We  use  these  facts  in  the  following  manner.  First  consider 
that  we  have  a fixed  coordinate  system  in  which  r and  r'  are 
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radius  vectors  to  the  observation  point  and  source  point 
respectively.  Next,  consider  that  we  have  another  coordinate 
system  in  which  the  radius  vector  to  the  observation  point 
is  denoted  as  r^  and  the  radius  vector  to  the  source  point  is 
denoted  r ^ . The  subscript  D is  employed  to  indicate  that 
the  second  coordinate  system  is  oriented  in  a special  manner 
with  respect  to  the  delta  function  that  appears  in  equation  49. 
The  relationship  between  these  two  coordinate  systems  will  now 
be  described  with  the  aid  of  the  following  figure. 


Figure  2.  Coordinate  Systems 

The  origin  for  both  coordinate  systems  is  on  the  interface 
and  z and  z^  are  measured  positive  as  the  distance  normal  to 
interface  increases  into  the  lossless  medium.  The  reason 
that  the  subscript  D is  employed  can  now  be  explicitly 
stated  as 


-D 


(50) 


i.e.,  D refers  to  the  fact  that  in  this  coordinate  system, 
the  origin  lies  directly  under  the  delta  function  source. 
The  reason  for  introducing  the  D coordinate  system  is  that 


r’jRT 
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and  satisfy  the  described  requirements  for  equations 
33  and  41  to  be  a valid  representation  in  that  system.  Upon 
using  this  special  representation  for  r^  we  will  find  that 
the  infinite  sums  in  equations  33  and  41  reduce  to  finite 
sums.  After  obtaining  the  simplified  representations  for 

(r^ , r^)  and  Ks (r^ , r^)  we  are  not  in  a position  to  use 
the  representations  in  equations  45  and  47  because  r^  is  not 
free  to  vary  in  the  D coordinate  system.  We  obtain  our 
simplified  representations  that  can  be  used  in  equations  45 
and  47  by  explicitly  performing  the  substitutions  implied  by 
the  following  equations 


G^^Cr.r’)  = G^i;)(rD(r,r' ) , r£(r,r'))  (51) 

Ks(r,r')  = Ks(rD(r,r’),  r£(r.r'))  (52) 


-4lA<'hl’-D') 


\ ih^z’ 


^olA^hl’-D^  2 e 1 ax 

^oA(hl’-D)  71  6 1 a2 


(56) 


(57) 


(58) 


h1A 

2ko 

ihiz1  - 
e ax 

(59) 

hxA 

ihiz’  " 
e L a 

(60) 

2ko 

y 

^ol  A (hl  ’ — 1 


Accordingly  equations  33  and  41  reduce  to  the  following 
finite  sums  in  the  D coordinate  system 


0°  / 

4^  y 2|a(^elA(hl’-D)  ^elA(hl’-D'> 


+ ^olA(hl’-D)  -olA^l’-D^ 

+ b(l  Whl^D>  KoX^l'-V 

+ NelA(hl’-D)  NelA(hl'^D) 

+ ^olA^l’^D^  1 A 1 ’ — D ^ ) 


(61) 
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MqIx^i-Ed)  p 


Jl(XpD)_  „ * iMz  2 

COS  4>j^  6 apjj 


+ - ^V^d’)  si"  ♦d  'lhlZ  SD 


(65) 


WVed>  - t Jc(iV  elhiz  “*  - 1 r-  V^d’  •~1~  SP 


hxX 


ih,  z ~ 
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— e 1 A 1 ’ — D ^ " t~  Jl(ApD>  cos  elhlZ  az 
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^1A(W  * r"  Jl(XpD>  sin  *0  elhlZ  az 
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^i(xPrj^  ihiz  - 
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To  obtain  G,  (r,r’)  and  K (r,r')  we  use  prescriptions  51  and  52 
J s s 

in  conjunction  with  transformation  equations  53  through  55  and 


Pn  = 


Pn  = 


IP  * P 

P - P' 


D P 
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az  * ap. 


(69) 
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We  are  now  in  a position  to  rewrite  the  EFIDE  and  MFIE 
operators  as  the  sum  of  three  terms  as  we  stated  in  the  abstract. 

The  first  term  describes  free-fieid  interaction,  the  sum  of  the 
first  two  describes  interaction  above  a perfectly  conducting 
ground  and,  naturally,  the  sum  of  all  three  terms  describes 
interaction  above  a finitely  conducting  ground.  The  corresponding 
decomposition  for  the  source  term  will  be  given  in  the  next  section. 

Equations  31  and  37  represent  the  first  stage  of  the  decomposi- 
tion, i.e.,  the  free-field  interaction  term  has  been  separated 
from  the  scattering  form.  Thus  we  must  further  decompose  the 
scattering  term.  Recalling  equations  61,  62  and  subsequent 
equations  63  through  68  we  can  cast  G (r^ , r^)  and  K^r^r^)  into 
the  following  forms 


G,  (rn,  rl ) — G,  a a + G . a . a,  +G-.  a a 
3s  -D'-D'  3sop  pD  p^  3sp4>  4>D  3spz  pn  : 


+ G,  a a + G,  a a 
3szp  z pD  3szz  z z 


( to  : 


Ks(-D'^D} 
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D D 
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+ K a a 

spz  pd  z 


( 7i : 


where 
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(72h) 


(73a) 


(73b) 
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The  above  integrals  are  well  studied  and  can  be  found  in  reference 
3 where  the  same  symbols  are  used  with  a slightly  different  but 
equivalent  representation  for  the  integrands. 
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To  proceed  with  the  decomposition  we  consider  the  limit  as 
conductivity  a of  the  ground  becomes  infinite.  One  can  show  that 


lim  V 
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and  equations  72  in  this  limit  becomes 


where  the  superscript  P means  perfectly  conducting.  By 
regrouping  terms  in  equations  72  and  utilizing  equations  74, 
we  obtained  the  following  decomposition 


«p  =p 

where  the  elements  of  G7  and  K are  given  by  equations  74 

=*=P  3 s ^ p s 

and  the  elements  of  G..  and  K (F  stands  for  finitely  conducting) 

j s s 

are  given  by 
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Thus  by  invoking  equations  31,  37,  75  and  76  we  can  complete  the 
promised  decomposition: 


- =-  =.p  —F 

G,_  = G + G,  + G, 
3s  o 3s  3s 


(78) 


**  — =p  =F 

K=  K+  + K* 
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where  Gq,  Kq  are  given  by  equations  32  and  38  respectively,  the 
elements  of  GP  and  KP  by  equations  74  and  the  elements  of 

=~F  =p  ■is  s 

G^s  and  Kg  by  equations  77. 

As  a check  to  the  above  calculations  we  should  arrive  at 
equations  74  by  employing  the  standard  forms  for  GP  and  KP. 
These  represent  the  scattered  parts  for  interaction  above  a 
perfectly  conducting  ground  and  are  given  by  the  well  known 
expressions 


-z')  I + 2S  a g (z* 
“ z z o 


(80) 


(81) 


where  gQ  is  the  free  space  scalar  Green's  function  given 
by  equation  33,  G2^  is  defined  by  equation  73a  and  Ir  is  a reflec- 
tion operator  given  by 


a 

z 


It  can  be  shown  by  some  rather  lengthy  algebraic  manipulations 
that  equation  80  is  identical  to  equations  74a  through  74d  and 
equation  81  is  identical  to  equations  74e  and  74f. 


SECTION  V 


REPRESENTATIONS  FOR  E (r)  AND  H,p(r) 


Using  equation  16  in  conjunction  with  the  coordinate 
notation  changes  described  to  obtain  equation  27  and  equation 
29,  we  obtain 


ET(r) 


iu)U0 


(r , r ' ) 


J(r')  dV’ 


(82) 


and  using  equations  35,  36,  and  82  we  obtain 


HT(r) 


K (r , r 1 ) 


J ( r ' ) dV ’ . 


(83) 


As  discussed  in  the  previous  section,  ET  (r ) and  H_T(r)  have 
meaning  in  their  own  right,  independent  of  their  role  in  the 
appropriate  EFIDE  or  MFIE . They  are  the  total  electric  and 
magnetic  fields  at  any  point  r in  the  lossless  medium  due  to 
a specified  source  distribution  J(r'),  r'  € and  equations 

82  and  83  include  all  interaction  with  the  lossy  half  space. 

The  work  presented  in  the  previous  section  enhances  the  utility 
of  equations  82  and  83  in  that  the  form  to  be  used  for  (r,r‘) 

is  the  simplified  form  given  by  equation  78,  while  the  form  to 

be  used  for  K(r,r')  is  the  simplified  form  given  by  equation  79. 

If  we  now  focus  our  attention  on  the  role  of  ET(r)  and  H^r) 
in  the  EFIDE  and  MFIE  then  we  consider  two  cases.  The  first  case 
is  when  we  desire  to  calculate  the  surface  current  density 
induced  on  the  scattering  body  and  the  source  of  excitation  is 
a specified  distribution  of  current  that  exists  somewhere  in  the 
lossless  medium.  For  this  case  we  must  employ  equations  82  and 

83  explicitly  in  conjunction  with  the  EFIDE  and  MFIE.  The  second 


case  is  one  of  considerable  interest.  This  is  the  case  when  the 
source  of  excitation  is  an  incident  plane  wave.  For  this  case 
we  make  use  of  the  fact  that  ET (r ) and  H_, (r ) are  simply  the 
incident  plus  reflected  fields  from  the  lossy  half-space  and 
we  do  not  have  to  explicitly  evaluate  equations  82  and  83. 

For  completeness  we  will  present  explicit  representations  for 
this  case.  These  representations  are  obtained  in  the  standard 
manner  in  which  the  incident  field  is  decomposed  so  that  it 
is  the  superposition  of  two  component  fields.  One  component 
field  has  its  electric  field  polarized  perpendicular  to  the 
plane  of  incident,  and  the  other  component  field  has  its 
magnetic  field  polarized  perpendicular  to  the  plane  of  inci- 
dence. The  analysis  leading  to  the  results  is  straightforward 
after  the  described  decomposition  is  employed.  For  an  inci- 
dent field  given  by 


E . = e E eik°n°*E 
— 1 o 

(84) 

with  e 

being  a unit  vector  along  the  polarization  directi 

on  and 

nn  being  a unit  vector 

in  the  propagation  direction,  the 

total 

fields 

are 

ET  = Ei  + E?  + EF 

(85) 

and 

+ + 

(86) 

where 

H.  = (E  /Z  ) 
-i  o'  o 

(nQ  x e)  ei]WE 

(87) 

f - -E0[<* 

• h)  h + (e  • v)  (n_  * h)  . e''"k°n^  - 

K J 

(88a) 
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EF  = E 


(Rjj  + 1)  (e  • h)  h -(R^  - 1)  (e  • v)  (nR  x h)j 
ik0n^*r  (88b) 


H = -(E  /Z  ) (e  • h)  (n  x h)  - (e  • v)  h 

— O O K 


cl  ik^rwr 

hje  o R - 


(88c) 


Hp  = (Eq/Zq)  [(Rh  + 1)  (e  • h) (nR  x h)  + (Rv  - 1) (e  • v)h] 

elkonR,£  (88d) 
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R and  R,.  are  the  usual  Fresnel  reflection  coefficients 
H v 

so  the  square  root  is  defined  to  have  a positive  real  part. 

We  note  that  as  the  conductivity  of  the  half-space  approaches 
infinity,  then  R -*■  -1  and  R,.  - 1.  It  can  now  be  seen  that 
in  this  limit  Er  and  vanish  and  this  was  the  reason  for  the 
decomposition  of  and  as  expressed  in  equations  85  and  86. 
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The  first  two  terms  in  these  expressions  represent  the  total 

fields,  incident  plane  reflected  due  to  the  presence  of  a 

F F 

perfectly  conducting  half-space.  The  addition  of  E and  H 

to  the  "perfectly  conducting"  total  fields,  yields  an  exact 

T T 

representation  for  the  total  fields,  E and  H , resulting  from 
the  presence  of  the  finitely  conducting  half-space. 
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SECTION  IV 


COMPARISON  BETWEEN  NUMERICAL  AND  EXPERIMENTAL 
DATA  FOR  METALLIC  OBJECTS  IN  FREE  SPACE  AND 
ABOVE  A PERFECTLY  CONDUCTING  GROUND 


In  this  section  we  compare  numerical  data  derived  via  our 
Magnetic  Field  Integral  Equation  (MFIE)  Code  to  experimental 
data  obtained  at  the  University  of  Michigan.  The  quantity  of 
interest  is  the  magnitude  of  the  current  density  induced  on 
metallic  objects  in  free  space  and  in  the  presence  of  a per- 
fectly conducting  ground  by  a monochromatic  plane  wave.  In 
particular  we  considered  a perfectly  conducting  circular  cylin- 
der in  free  space  and  also  above  a perfectly  conducting  ground 
and  a perfectly  conducting  aircraft  model  in  free  space.  Thus 
figure  3 presents  the  comparison  for  the  magnitude  of  the 
induced  current  density  as  a function  of  kh  (k  being  the  free 
space  wave  number  and  h the  halflength  of  the  cylinder)  at  the 
two  points  P and  Q indicated  on  the  graphs.  The  incident  wave 
is  polarized  with  its  electric  vector  parallel  to  the  axis 
and  the  propagation  vector  k is  perpendicular  to  the  axis. 

The  agreement  is  generally  very  good.  Notice  that  the  peak 
value  obtained  via  the  MFIE  code  is  higher  than  the  one 
measured  at  the  University  of  Michigan.  This  is  in  agreement 
with  the  results  obtained  by  Sassman  (ref.  9)  who  calculated 
the  total  current  I as  a function  of  kh.  Sassman  obtained 
84a  for  the  peak  value  of  the  I/HQ  and  this  should  be  smaller 
than  27ra|Jp|/Ho  since  |jp|  is  maximum  at  P.  The  code  gives 
2iTa|Jp|%89a  and  the  measured  data  gives  84a.  (Figures  4 and  5 
are  taken  from  the  University  of  Michigan  report  since  the 
solid  curves  in  figure  3 are  our  smoothed  out  drawings  for 
the  original  graphs  in  figures  4 and  5 and  as  such  they  are 
subjective . ) 
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Figure  3a, b.  Comparison  Between  "Free  Cylinder"  Experimental  Results 
at  University  of  Michigan  (Solid  Curves)  and  TDK's 
Code  ( Dots) . 


In  Figure  6 the  cylinder  is  situated  above  a perfectly  conduct- 
ing ground  with  its  axis  parallel  to  the  ground.  The  incident 
field  is  the  same  as  in  Figure  3 and  the  distance  of  the  cylin- 
der axis  from  the  plane  is  5a. 

Again,  Figure  7 refers  to  the  original  graph  in  the  Univer- 
sity of  Michigan  report.  Figure  8 depicts  the  same  situation  as 
in  Figure  6 except  that  the  cylinder  axis  is  much  closer  to  the 
ground.  (Figure  9 is  the  corresponding  graph  in  the  University 
of  Michigan  report.)  Notice  that  the  deviation  near  the  two 
maxima  is  much  greater  than  everywhere  else.  To  satisfy  our- 
selves that  the  discrepancy  between  the  experimental  values 
of  current  density  and  the  corresponding  numerical  predictions 
were  not  the  fault  of  the  finite  number  of  zones  used  to  parti- 
tion the  surface  of  the  cylinder,  we  systematically  increased  the 
number  of  zones  until  we  obtained  what  appeared  to  be  a con- 
verging sequence.  For  the  cylinder  in  free  space  and  for  the 
cylinder  5 radii  above  the  perfectly  conducting  ground  plane 
we  increased  the  zoning  until  for  each  frequency  the  response 
was  insensitive  to  zoning  changes.  However,  the  small  change 
in  the  predicted  resonant  frequency  due  to  a change  in  zoning, 
coupled  with  the  sharp  resonance  for  the  case  in  which  the 
axis  of  the  cylinder  was  1.5  radii  above  the  ground  plane, 
caused  the  above  convergence  criterion  to  be  too  stringent  to 
be  practical  for  frequencies  near  this  resonance;  instead,  in 
this  region,  we  required  that  both  the  resonant  frequency  and 
the  peak  value  be  insensitive  to  zoning  changes. 

The  zoning  scheme  that  produced  the  results  found  in 
Figure  8 was  18  zones  along  the  length  and  10  zones  around 
the  circumference,  while  for  Fiqure  6 it  was  sufficient  to  use 
10  zones  along  the  length  and  6 zones  around  the  circumference. 

In  bott  cases  there  was  one  zone  along  the  radius  of  the  endcap. 


Figure  6.  Comparison  Between  "Cylinder  Parallel  to  a Perfectly  Conducting 
Ground"  Experimental  Results  at  University  of  Michigan  (Solid 
Curve)  and  TDR ' s Code  (Dots  for  d = 5a. 
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Figure  9.  The  University  of  Michigan  Original  Plot 
for  the  Solid  Curve  in  Figure  8. 


Figure  10  is  a schematic  for  the  aircraft  model  considered 
at  the  University  of  Michigan  the  exact  dimensions  of  which  are 
given  in  Table  1.  The  experiments  were  conducted  with  the  brass 
scale  model;  however,  our  graphs  refer  to  the  unsealed  aircraft 
dimensions  given  in  Table  1.  Figures  12,  14  and  16  present  the 
comparison  for  the  magnitude  of  the  induced  current  density 
along  the  top  of  the  fuselage  for  kh  = 0.82,  1.694,  20.14,  where 
h is  the  fuselage  halflength.  Figures  13,  15,  and  17  refer  to 
the  bottom  of  the  fuselage.  The  incident  plane  wave  has  its 
electric  vector  polarized  parallel  to  the  fusalage  and  the 
propagation  vector  perpendicular  to  the  fuselage  axis.  The 
distance  £(in  meters)  shown  in  the  graphs  is  the  arclength 
defined  in  Figure  11.  Notice  that  the  deviation  between  experi- 
mental and  theorectical  values  on  the  fin  (vertical  stabilizer) 
is  much  greater  than  on  the  fuselage.  (In  figure  13  our  results 
were  off  scale  and  are  not  shown.)  A possible  explanation  for 
this  discrepancy  may  be  based  on  the  fact  that  the  measuring 
sensors  are  of  finite  size  and  do  not  pick  up  the  fields  right 
on  the  surface.  This  is  important  because  the  elliptical  cross 
section  of  the  fin  has  a ratio  av/bv  % 6 , i.e.,  the  curvature 
about  the  points  where  the  measurements  were  performed  varies 
rapidly  and  consequently  the  fields  can  change  rapidly  as  we 
move  away  from  the  surface.  This  is  at  least  true  for  the 
special  case  of  a perfectly  conducting  cylinder  of  infinite 
length  and  elliptical  cross  section  immersed  in  a magnetostatic 
field  parallel  to  the  minor  axis.  For  this  case  we  found  that 
at  a distance  (in  the  direction  of  the  major  axis)  equal  to 
0.25  bv^  0.025  in  (bv  is  the  minor  semi-axis)  the  azimuthal  com- 
ponent of  the  total  magnetic  field  is  equal  to  half  its  value 
on  the  surface.  The  above  discussion  indicates  that  our 
higher  values  on  the  surface  of  the  fin  may  actually  be  compat- 
ible with  the  lower  values  obtained  with  the  finite  size  sensors 
which  of  necessity  measure  fields  at  points  off  the  surface. 


Table  1. 

MODEL  SPECIFICATION 


Figure  12.  Bottom  of  Fuselage,  kh  = 0.82.  (h:  fuselage  half 

length) . 


Figure  13.  Top  of  Fuselage,  kh  = 0.82.  (h:  Fuselage  half  length. 


Figure  15.  Top  of  Fuselage,  kh  = 1.694.  (h:  Fuselage  half  length. 


Figure  16.  Bottom  of  Fuselage,  kh  = 2.014.  (h:  Fuselage  half  length) 
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